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Abstract. We derive some new relationships between matrix models of Chern-Simons gauge 
theory and of two-dimensional Yang-Mills theory. We show that g-integration of the Stieltjes- 
Wigert matrix model is the discrete matrix model that describes g-deformed Yang-Mills theory 
on S^. We demonstrate that the semiclassical limit of the Chern-Simons matrix model is equiv- 
alent to the Cross- Witten model in the weak coupling phase. We study the strong coupling 
limit of the unitary Chern-Simons matrix model and show that it too induces the Gross- Witten 
model, but as a first order deformation of Dyson's circular ensemble. We show that the Suther- 
land model is intimately related to Chern-Simons gauge theory on S^, and hence to g-deformed 
Yang-Mills theory on . In particular, the ground state wavefunction of the Sutherland model 
in its classical equilibrium configuration describes the Chern-Simons free energy. The corre- 
spondence is extended to Wilson line observables and to arbitrary simply-laced gauge groups. 



1. Introduction and summary of results 

Matrix models have been a subject of much interest in gauge theory for over three decades [1]. 
In this paper we study random matrix models [2] related to Chern-Simons theory [3], and both 
ordinary two-dimensional Yang-Mills theory 0] and its g-deformed counterpart [5l [6] . We will 
develop some new relationships between the Chern-Simons matrix models and some of the matrix 
models that appear in two-dimensional Yang-Mills theory. 

Let us begin by giving the definition of the joint probability distribution P{M) for the matrix 
elements of an x matrix M [2], 

(1.1) P{M) = Z-^^ exp ( - Tr V{M)) , 

with any potential V{M) such that the partition function exists. The integration of (jl.ip 
over parameters related to the eigenvectors of M leads to the well-known joint probability 
distribution of the eigenvalues [2] 

TV 

(1.2) p(xi,...,x^) =z^^ n n^^p(-^(^^)) • 

i<j 1=1 



The level repulsion in (ll.2p . described by the Vandermonde determinant, originates as the Ja- 
cobian of the transformation when passing from integration over independent matrix elements 
to integration over the smaller space of N eigenvalues. The integer /3 = 1,2,4 describes the 
symmetry of the ensemble (orthogonal, unitary and symplectic, respectively). 
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These symmetries have been extended using the Cartan classification of symmetric spaces [7j 
(see also [8l|9l[T0]). For example, the possible Jacobians of the transformations to radial coor- 
dinates are given by [8] 

(1.3) jW(x) = n (^")"" ' 

aGA+ 

j(-)(x) = Yl (a-i sinh(x«))"'" , 
aeA+ 

j(+)(x) = J] {a-' sin(x„))™" , 

for the various types of symmetric spaces with zero, negative and positive constant curvature, 
respectively, and for an arbitrary non-zero constant a. The products are taken over all positive 
roots of the restricted root lattice, with the multiplicity of the root vector a E A"*", and 
Xa '■= {x, a) = Xi a* are canonical coordinates on a maximal abelian subalgebra of the 
tangent space. 

The case of the A^-i root system and zero curvature leads to the well-known Gaussian matrix 
model distribution (we write only the Hermitian case) 

N 

(1.4) P(^\xi, ...,xn) = Z^'1[ {x, - Xjf llexp{- xj/2) , 

i<j i=l 

whereas the same root lattice but in the case of negative curvature yields 

N 

(1.5) P^-\xi,...,XN) = Z^'llsmh^a{xi-Xj)) J] exp ( - xf /2) . 

i<j 1=1 

In the B-model topological string theory, this is the matrix model that describes U{N) Chern- 
Simons gauge theory on [12] if a = y/g^, where gs is the string coupling constant which 

is related to the usual integer Chern-Simons level k by 

(1.6) gs = . 

^ ' ^ k + N 

The exponential mapping Xi = logUj/a brings (|1.5|) into the usual form p.2|) (with (3 = 2), and 
hence one can apply the orthogonal polynomial method of random matrix theory [2] to solve the 
model |13[ll4j. The orthogonal polynomials here are the Stieltjes-Wigert polynomials |13 [ I14 [ [T5] 
or, if we work with a unitary matrix model [TH], the Rogers-Szego polynomials |15| . 

Chern-Simons theory on is equivalent to g-deformed Yang-Mills theory on 5^, with the 
identification q := e~^^. This result can be derived at weak coupling either directly via lo- 
calization of the three-dimensional gauge theory by regarding as a Seifert manifold through 
the Hopf fibration — ?■ S'^ \17\ I18j . or by recasting the two-dimensional gauge theory as a 
sum over instantons and explicitly demonstrating its equivalence with the Chern-Simons matrix 
model [ini [20] . At weak coupling, the two-dimensional U (N) gauge theory thus reproduces the 
perturbative A-model topological string partition function Ztop for the resolved conifold geom- 
etry in the large N limit. The Stieltjes-Wigert polynomial is the average of the characteristic 
polynomial in the matrix model, and hence describes B-brane amplitudes on the conifold [21J. 

In the next section we show that (^-integration of the Stieltjes-Wigert matrix model directly 
gives the discrete matrix model that describes the strong coupling expansion of g-deformed two- 
dimensional Yang-Mills theory. In the large limit, the strong coupling series has zero radius 
of convergence |19) . This is in sharp contrast to ordinary (undeformed) Yang-Mills theory on 
S'^, which undergoes a third order phase transition at large A^ [22] and possesses a double 
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scaling limit which lies in the universality class of the Gross- Witten unitary matrix model 
The double scaling limit of the Gross- Witten model is also of interest in the study of unitary 
matrix models of string theory in zero dimensions [23], of the solution of SU{2) Seiberg- Witten 
theory [25], and of Type OA and OB string theories in one dimension ^B] where it is argued to 
describe the universality class of pure two-dimensional supergravity. 

In the next section we will demonstrate that both the weak and strong coupling limits of the 
finite Chern-Simons matrix models are also governed by the Gross- Witten model [27]. Recall 
that this is the unitary one-matrix model which arises as the one-plaquette reduction of the 
combinatorial quantization of Yang-Mills theory in infinite spacetime. In two dimensions the 
reduction is exact and described by the partition function [J^ 



(1.7) Z^^{a):= f dC/exp(-aTr(C/ + [/^ 

Ju(N) ^ 

•^0 i=l i<j ^ 

where dU denotes the bi-invariant Haar measure for integration over the unitary group U{N). 
This matrix model belongs to the class of symmetric spaces associated with the A^v-i root 
system and positive curvature. 

We will thus show that the g — )• 1 limit of the unitary Chern-Simons matrix models are 
equivalent to two-dimensional Yang-Mills theory on the plane M?. Identifying a = = 
and using the exact integration of (|1.7|) in the weak-coupling phase |28t eq. (5.25)], we can write 
down a simple all-orders expression for the U{N) Chern-Simons partition function Zj\[^k{S^) on 
in the weak-coupling limit as 

nI,N, ) (:.(-l)-/=e-(-./..) . 

This exact expression should reproduce the topological string partition function Ztop in the large 
N limit to leading orders of perturbation theory. 

As the unitary one-matrix model (jl.7p undergoes a phase transition to a strong coupling 
phase in the large N limit, while the q-deformed gauge theory on S"^ always remains in its weak 
coupling phase, the nature of this relationship must be drastically different at strong coupling. 
We will show that the g — )• limit of the unitary Chern-Simons matrix model is equivalent 
to Dyson's circular ensemble, together with an infinite tower of higher Casimir deformations 
corresponding to multicritical extensions |24J, the lowest order of which is described by the 
Gross- Witten model. These deformations are analogous to those which arise in Yang-Mills 
theory on the noncommutative torus \29l |30j . This feature puts the g-deformed gauge theory 
into the context of noncommutative deformations of Yang- Mills theory [31] . 

In the last section we show that the celebrated Sutherland model [32] . a central model in the 
theory of one-dimensional integrable systems, is also directly related to Chern-Simons gauge the- 
ory. This provides another connection between two-dimensional Yang-Mills theory and Chern- 
Simons theory, and moreover between integrable models and Chern-Simons theory on the three- 
sphere S^. The Sutherland model is the exactly solvable system on a circle of circumference L 
defined by the A-particle quantum Hamiltonian operator 



^ 52 2A(A- l)7r2 ^ / 7r(g, -g,) 
I ^ ' sm 



i=l * i<j 
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with the ground state wavefunction 




) 



A 



(1.10) 



i<j 



The relationship foUows immediately when the model is considered in a simple fixed crystalline 
configuration. As we will discuss, this is not an artificial configuration, and it corresponds to 
the classical equilibrium state of the Sutherland model which is naturally reached in the strong 
coupling limit A — )■ oo. We will show that the probability density distribution, evaluated in an 
equally spaced configuration, reproduces the partition function for Chern-Simons gauge theory 
on with gauge group SU{N), which is one of the simplest quantum topological invariants 
that can be obtained from Chern-Simons theory [3]. We demonstrate that this correspon- 
dence extends to Wilson line observables, and hence to other quantum topological invariants, 
by considering excited states of the Sutherland model associated to particular non-equilibrium 
configurations of the particles. It also holds between generalized Calogero-Sutherland models 
and Chern-Simons theory with other simply-laced gauge groups (associated to orthogonal and 
symplectic ensembles). 

Recall that the Chern-Simons partition function Zjv,fc(5''^) (or equivalently the free energy 
Fjv,fc(S'^) = log Z7v,A:(5'^)) is central in the study of topological strings, for example. The matrix 
model description is useful in topological string theory, and its physical applications have been 
exploited mainly in that context thus far. Although there have been extensive applications of 
Chern-Simons theory to condensed matter physics (most notably to the fractional quantum Hall 
effect), our application here deals with a many-body description of Chern-Simons theory in its 
non-abelian version, at least for certain simple three-manifolds such as S^. Given the natural 
appearance of random matrix theory, and its connection with exactly solvable models in one 
dimension [Rj and with Laughlin wavefunctions [33] , it seems worthwhile to further explore this 
many-body description in itself and its possible role in condensed mater physics. 

This correspondence also gives an intriguing new relationship between two-dimensional Yang- 
Mills theory and certain one-dimensional integrable systems, which is rather different in spirit 
from previous relationships [341 [35j that considered wavefunctions of Calogero-Sutherland and 
Calogero-Moser models as reductions of those for the undeformed gauge theory on a cylinder 
with appropriate Wilson line insertions. The g-deformed gauge theories on S'^ may thus provide 
computationally useful means for exploring various aspects of one-dimensional exactly solvable 
models. In particular, as the large deformed gauge theory on exhibits no phase transitions, 
there is a natural and simple well-defined large N field theory limit of the Sutherland system 
which is equivalent to a generalized Calogero-Sutherland model based on the su{oo) Lie algebra. 
In [31] this limit is instead realized by replacing SU{oo) with the centrally extended SU {N) loop 
group, which is related to Yang-Mills theory on an infinite cylinder with Wilson line insertions. It 
would be interesting to better understand the explicit relationship between the two approaches, 
given the known equivalence between the partition functions of q-deformed Yang-Mills theory 
on 5^ and of ordinary Yang-Mills theory on a cylinder with trivial holonomies around the two 
boundary circles of the cylinder [36] . 

Acknowledgments. This work was supported in part by grant ST/G000514/1 "String Theory 
Scotland" from the UK Science and Technology Facilities Council. MT thanks Mark Adler 
for hospitality at the Mathematics Department at Brandeis University, and the Department of 
Mathematics at Heriot-Watt University for a productive stay. 
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2. From Chern-Simons theory to the Gross-Witten model 

2.1. Another derivation of g-deformed Yang-Mills theory. Let us begin by computing 
the g-integration of the Stieltjes-Wigert matrix model. We will use Jackson's integral, which in 
the single variable case is given by [371 EB] 

/"OO ^ 

(2.1) / d,uw{u) = {l-q) ^ wiq^q^. 

•'O n=-oo 

If the function w{u) is the weight function of the Stieltjes-Wigert matrix model, i.e. the log- 
normal distribution w{u) = exp (— log^ u/2gs), then the g- integration gives 



(2.2) 



poo 

/ dqu e-i°g'"/29. = J2 , 



with the usual identification q = e The right-hand side of (j2.2p is proportional to the 
theta-function 



(2.3) eoo(^k)= X; q 



at the particular value z = q. This theta-function is also the weight function of the unitary 
matrix model that describes Chern-Simons gauge theory [16\ 115] , We shall see below how this 
feature relates the Chern-Simons matrix model with the Gross-Witten model. 

Let us now compute the g-integration of the log-normal weight function in the multivariable 
case. For the Stieltjes-Wigert matrix model, we then have 



(2.4) 1^ ^-log'u,/2gs Yl (, 



\Ui - u 



i=l i<j 

N 



which coincides with the partition function of the g-deformed two-dimensional Yang-Mills model 
on 5^ with gauge group U{N) [5l |6]. Recall that, using solutions of the log-normal moment 
problem, the equivalence between the continuum and discrete matrix models was demonstrated 
in this case in ^\ , and it reads 



(l)""7:nte-''-n(-.n.(2^)) 

1=1 i<j ^ ^ 

= ^;v(^) e-tE.n? -Q (2sinh(|(n,-n,))) 



where 

(2.6) AM ■-- 



r!,i,...,njv=— oo i<j 



^_(l_2Ar+3Ar2)/2 \^ 
(- 9^/^-^)00 (-^^-^/^)oo('?)oo J 



with (g)o := 1, {q)m ■= ni<j<m (l " Q^) for rn £ N, and {q)oo := HjeN (l " Hence 
(/-integration and ordinary integration of the Stieltjes-Wigert distribution gives the same result 
(up to overall normalization), the U{N) Chern-Simons partition function on S^. 
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2.2. Weak-coupling limit. We will now demonstrate that the Chern-Simons matrix model 
yields, in the weak-coupling limit (7^ — t- (g = e"^" — )■ 1), the Gross- Witten model ()1.7p with 
coupling g = — )• 0. For this, we note that the distinctive log-square behaviour of the weight 
function w{u) = exp iy—\o^ u/2gs) for the Stieltjes-Wigert matrix model also constitutes the 
natural definition of g-growth [4Ql [H] . We use results of ^40j to study its relationship with the 
different definitions of g-exponential functions, and in turn to better understand its g — )• 1 limit. 
The crux of our analysis is that the weight function of the Stieltes-Wigert matrix model can be 
given in terms of g-exponential functions, but without the usual rescaling u — )■ {1 — q)u that is 
necessary to recover ordinary calculus from g-deformed calculus when g — )• 1 |371 [38] . 

As we mentioned in the previous section, the matrix models that appear in Chern-Simons 
theory can be solved exactly by using orthogonal polynomials. These orthogonal polynomials 
are completely characterized by a discrete scaling symmetry satisfied by the weight function [39] 

(2.7) w{qu) = ^/quw{u) . 

This self-similarity property is called the g-Pearson equation, and it uniquely determines the 

orthogonal polynomials up to normalization |39j . 

Let us begin by briefly reviewing g-deformed expressions for exponential functions and how 

they tend to their classical counterparts in the limit 9 — )• 1 |37t 138] . For n G C, the two main 

definitions of a g-exponential function in the literature are given by 

00 j 

w 1 



(2-8) e,{u) = Y.T7^ 



Egiu) = 2^ 7-T = i-u; q)^ = (eg(-n)) , 

where (u; q)oo '■= YijeNo (l ~ ^"') • The connection with the usual exponential function is given 
by the functional equation 

(2.9) eq{u) Eq{-u) = 1 
and the weak-coupling limits 

(2.10) lim eg( (1 - q) u) = lim Eg{{l - q) u) = e"" . 

We now recall that the functional equation ()2.7p solves the matrix model [39]. Furthermore, 
one can show [IQ] that the entire function f{u) = Eq{u) satisfies (1 + u) f{qu) = f{u) and 
the entire function g(u) = eg-i{—u~'^) satisfies ug{u) = {1 + u) g{u). It follows that the entire 
function 

(2.11) h{u) = f{u)g{u) = {-u;q)^ (^-qn-^;q)^ 

is a solution of the g-difference equation uh{qu) = h{u). By replacing q — )• q~^ and rescaling 
u ^ qu, the functional equation for h{u) reads 

(2.12) h{qu) = quh{u) . 

This is essentially the g- Pearson equation ()2.7p . but with q instead of y/q appearing the right- 
hand side, and hence numerical prefactors will be different if we use (j2.12p instead of (j2.7p as 
the defining characterization of the Stieltjes-Wigert weight function w{u). 
Thus the explicit g-deformed expression that satisfies ()2.12p is given by 

(2.13) w{u)=Eg-i{qu)eq{-{quy^) . 
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An elementary manipulation shows Eq-i{u) = eg{u), and so from (|2.1Up the weak-coupling limit 
(7 —7- 1 is given by 

(2.14) lim w(u) = exp ( | exp ( | = exp ( (u + u^^] 

In the unitary Chern-Simons matrix model, we substitute u = e'^ with < < 27r to get 

(2.15) lim w(u = e = exp f cos^ | , 

and we arrive at the Gross- Witten model ()1.7p with a oc j^. Since g — )• 1 (g^ — ^ 0) we can 
write a ~ J- . The factor of 2 in (|2.15p is not accurate since we have used (|2.12p instead of (|2.7p . 

The same result is contained in the g — >• 1 limit of the (^-Hermite polynomials studied by 
Ismail, Stanton and Viennot in [l2]. They consider polynomials with the same orthogonality 
properties as the ones that solve the Chern-Simons matrix model (see the discussion of the g — t- 
limit below), but then apply a rescaling of the variable x — t- — g/2) x in order to obtain 
the correct q 1 limit in their weight function v{x,q) as mentioned above. If we undo this 
rescaling and take into account that with their definitions x = cos (0/2), then their result 

(2.16) lim v(x, q) = exp i-x^ 12) 

leads in our case to 

^■ / \ 2 o6\ _i/n--„-) / COS0 

(2.17) lim v(x,q) = exp cos - = e ^' exp 

q^i V 1-9 2/ V 1-9 

2.3. Strong-coupling limit. Let us now study the connection between the Chern-Simons ma- 
trix model and the Gross- Witten model in the strong-coupling limit ^ oo {q = e~^^ — )• 0). 
Recall first that both the Stieltjes-Wigert polynomials, defined on the positive real half-line, 
and the Rogers-Szego polynomials, defined on the unit circle |43| . solve the matrix models that 
appear in Chern-Simons theory \14\ [T5] . In the Hermitian case this was demonstrated in [T3] , 
while the unitary matrix model was proposed in [16] (but without consideration of the associated 
Rogers-Szego polynomials). The explicit connection between the two models was demonstrated 
in [15j. The direct relationship between the two systems of orthogonal polynomials is well- 
known [Uj. Therefore, the partition function for Chern-Simons gauge theory on with gauge 
group U {N) can be represented both as the Hermitian matrix model 



(2.18) 



ZN,k{S^)= j [dM] exp(^-^Tr(logM)2^ , 



with [dM] the natural invariant measure for integration over the space of N x N Hermitian 
matrices, and as the unitary matrix model [T6] 

(2.19) ZN,k{S^)=[ dU deteoo{U\q) , 

JU{N) 

where the theta-functions 0oo(6 ^^'\q) are defined on eigenvalues of unitary matrices U by (j2.3p . 
The Rogers-Szego polynomials are defined by [43j 

(2.20) H^{z\q) := M , M := , , • 

fr'oUlg b\g (9)i(9)n-J 

They satisfy an orthogonality relation on the unit circle given by 

(2.21) / ,^M- 9"'/'^ \q)M- 1''^' ^ I 9) eoo(^|9) = % W . 
J\z\=i ^TTiz q 
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The orthogonality coefficients hm = are those of the Stieltjes-Wigert polynomials, that lead 
to the Chern-Simons partition function |14|i^ 

The weight function appearing in ()2.2ip is u}{z) = @QQ{z\q) which can be expanded in powers 
of q and z using ()2.3p . Since here z = e^^ lives on the unit circle, the expansion reads 

(2.22) a;(cos 0) = 1 + 2 ^ cos 6^ + 2 ^ g"'/^ cosn9 . 

n>2 

The first two terms in (j2.22p give the first order approximation of the weight function for the 
Chern-Simons matrix model as g — t- 0. This result is very different from the — >• 1 behaviour. 
The leading constant term in ()2.22p corresponds to Dyson's circular ensemble |2]. The first 
order correction for small q is given by the weight function of the Gross- Witten model (|1.7p with 
a = ^Jq = e~^°/^ — t- 0. Corrections to this term are given by higher powers of cos^. 

These behaviours are similar to those obtained in [35| for a unitary matrix model with a 
one-parameter family of weights which defines a g-deformation of the circular ensemble. The 
a = member of this family is the weight function [45\ eq. (3)] 

(2.23) vo{e,q) = \{^ e'';q)^f ■ 

In the limit g — )• it tends to Dyson's circular ensemble and in the limit g — )• 1 to the Gross- 
Witten model. This weight is very similar to that given by the theta-function above (although 
the two infinite product expansions are different). 

2.4. Comparison with multicritical and noncommutative deformations. If one keeps 
only a finite number of terms in the expansion (j2.22p . then the corresponding matrix model is 
the multicritical polynomial generalization of the Gross- Witten model discussed in [24] . These 
results have been generalized and unitary matrix models of the form 

„ oo 

(2.24) Z^"it(t) = / dC/exp(^ ^ tjTvW 

j=-oo 



U{N) 



are of interest in various areas of field theory (see for a recent review) . In general, the cou- 
pling parameters tj are restricted by solutions of Virasoro constraints on the partition function 
Z^"'*(t). However, if one chooses them to be the coefficients of a ^-series, and in particular 
those of the theta-function in ()2.22p . then one obtains the Chern-Simons matrix model. 

Our realization of the unitary matrix model describing Chern-Simons theory as a certain q- 
deformation of the Gross- Witten model is also reminiscent of what occurs in two-dimensional 
U{N) Yang-Mills theory on a noncommutative torus. In [29] it is shown that the weak-coupling 
limit of the partition function in this case can be regarded as coming from a modification of 
ordinary gauge theory by the addition of infinitely many higher Casimir operators to the action. 
The addition of higher Casimir operators to ordinary, two-dimensional Yang-Mills theory leads 
to generalized Yang-Mills theories [1]. 

The relationship between these noncommutative deformations and our g-deformations is most 
transparent in the combinatorial quantization of the gauge theory described in ^30j, which yields 
a discrete family of unitary matrix models parametrized by A^-th roots of unity C and partition 
functions given by [301 eq- (92)] 

(2.25) zJJ2(a) = V-i- / dU XRiU)exp(-aTr{CU + CU^ 

o dim-K Ju(N) ^ 



"'^ These are the same orthogonality coefficients that were considered in [15]. The g"™ term differs from the 
coefficients of [T3] , but it only contributes a phase to the partition function. See [TJ] and the discussion of framing 
in Section [3.21 b elow for more details. 
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where the sum runs over all irreducible unitary representations R of the gauge group U (N) 
with dimension dimi? and characters XR- The truncation of the sum in ()2.25p to the trivial 
representation is independent of C [281 §§5.2.1] and corresponds to the Gross- Witten reduction 
of the combinatorial quantization of ordinary gauge theory in two dimensions. The noncommu- 
tative gauge theory partition function thus generalizes that of ordinary Yang-Mills theory by 
perturbing it by a sum over non-trivial representations of the unitary group. Again, in contrast 
to the g-deformed gauge theory on S"^, this matrix model admits large N phase transitions 
and non-trivial double scaling limits which converge to the continuum (noncommutative) gauge 
theories [301 EI]. 

3. From the Sutherland model to Chern-Simons theory 

3.1. Derivation of the Chern-Simons free energy. Let us now establish the relationship 
between the Sutherland model and the Chern-Simons matrix model that we mentioned in the 
first section. We shall pay particular attention to the square of the ground state wavefunction 
(jl.lOp for the value A = ^ of the Sutherland coupling, which is given by 

Trfa -qj) 



(3.1) P{qi,...,qN;L) = {^o{qu---,qN;hW =11 



sm ■ 



L 

and evaluate this probability density in the fixed and equally-spaced configuration qi 
with c € M an arbitrary constant. Then we easily find 

vr (i - j) 



(3.2) P(c-l,...,c-iV;L) = J] 



sm ■ 



L 

t<j 

n/ . vr . 27r . tt{N yj . at. /vr . 

( sm- sm— • • • sm j = H sm [yJ 

This is proportional to the Chern-Simons partition function on [3j once we identify the 
arbitrary length parameter L as 

(3.3) L = k + N 

in terms of Chern-Simons parameters, the integer level k and the rank of the gauge group N. 
However, there is an additional factor (k + N)^^^"^ in front of the product term in the Chern- 
Simons partition function [3]. Thus we require a renormalization of the Sutherland wavefunction 
(jl.lOp by multiplying it with the parameter L~^/^ . Then one has 



(3.4) *o(9i,...,9^;iL):=^n (^"^^^^^ 



1/2 



with 

(3.5) {^o{l.....N-\,k + N)f = ZN,k{S^) ■ 

In the sequel we will drop the tilde notation on the renormalized wavefunctions. 

We could have obtained the same result directly from the ground state wavefunction (jl.lOp 
with the coupling A = 1. In this case the Sutherland model corresponds to a theory of free 
fermions on a circleU Furthermore, the configuration employed, with particles equally spaced 
on the circle, corresponds to the static classical equilibrium state of the Sutherland model. 
This condition is essential to preserve integrability when the model is discretized to construct 



connection between Chern-Simons theory and free fermions at finite temperature is discussed in [47) . 
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spin chain models (see [l^ for a review). It follows that we can also write the free energy of 
Chern-Simons theory on 5^ as 

(3.6) FN,k {S^) = lim \ log {qi, . . . ,qN; X,k + N) . 

A-s-oo A 

Recall that the static classical equilibrium configuration of the Sutherland model, together 
with an additional spin interaction term in the Hamiltonian, is the one that leads to the Haldane- 
Shastry spin chain model |49^ [50] 1^ Hence the part without the internal spin degrees of freedom 
naturally induces the Chern-Simons partition function. However, the translationally invariant 
spin wavefunction [511 eq. (5)] 

(3.7) ^o(.i,...,.^;«,iV)=5E„..,oe'tS«.o(«....-i) J| (sin^^^) 

l<m<n<N ^ ^ 

with s„ = lb 1 local spins and a a positive real number, can be expressed in terms of hard-core 
boson variables. When is even, one can identify the spin up (resp. down) states with empty 
(resp. occupied) hard-core boson states to write [5T| eq. (6)] 

(3.8) ^r(ni,...,n.;a,iV)=e-J^»"' n (^in^^i^^i^)'" , 

i<M<f 

ni>nj 

where Ui = 1, . . . ,N denote the positions of ^ hard-core bosons in the periodic chain. According 
to our analysis above, this also yields the Chern-Simons partition function on S^. 

3.2. Framing dependence. The many-body quantum system with ground state wavefunction 



(3.9) ^^^{q,,...,q^;9.,L) = ^lle-'^?/'^^ JJsinh^i^ 

i=l i<j 

is the one that appears in Chern-Simons theory [U], with the normalization constant giving 
the Chern-Simons partition function when L = 1 [T3]. It is also possible to obtain from ()3.9p the 
Chern-Simons partition function Z^^^iS^) in the same way as above. This was already carried 
out in [36j in the context of Brownian motion (see also [52\ Appendix B]). 
The relevant density considered in |36j is given by 

N 

(3.10) <°""(gi, . . . , czTv; t) = n n ' 

i=l i<j 

with t = —1/gs- The Gaussian factor in (j3.10p only contributes to a part of a framing factor [3l 
[53] when qi,. . . ,qN are taken to be the components of the Weyl vector p of U(N), which is 
given by 

N 



(3.11) p ^ 



E 

1=1 



N ^\ 



with Cj the standard basis of unit vectors in M^. This is simply the equilibrium configuration 
of the Sutherland model that we considered above. The contribution of the Gaussian factor is 
then the phase [3B] 



(3.12) &-\P\ /* = exp 



More generally, one can include internal SU{N) colour degrees of freedom for the particles. 
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Thus the model that leads directly to the Chern-Simons partition function with the canonical 
framing on is the Sutherland model, while the model with ground state wavefunction ()3.10p 
leads to the Chern-Simons partition function in the matrix model framing [363 . 

Generally, the contribution of a framing IT on a three-manifold M (here M = S^), i.e. a 
choice of trivialization of the bundle TM © TM, to the partition function for gauge group G 
(here G = U{N)) is parametrized by an integer s G Z and given by |53] 

^, /27ris \ /27risfcdimG\ /nislpl'^ k\ 

(3.13) J(M.n) = exp (— c) = e,p (— -j^) = .xp {j^) , 

where h is the dual Coxeter number of G (here h = N). In the second equality here we have 
used the explicit expression for the central charge c of the associated Wess-Zumino-Witten 
conformal field theory based on the affine extension of G. In the third equality we have used 
the Freudenthal-de Vries formula to relate the dimension of the Lie group G to the length of its 
Weyl vector p (see [53] for details). Thus if we write the framing factor (|3.13p in the form 

,3.M) ,(M.n).exp(^)exp(-^). 

we see that the contribution (j3.12p is given by the second factor in (j3.14p . 

A phase contribution that does not include the level k of the Chern-Simons gauge theory, 
as in the first factor of ()3.14p . appears when one goes from the hyperbolic to the trigonometric 
case and conversely. Indeed, one can also consider the hyperbolic Sutherland model with ground 
state wavefunction 

(3.15) <^P(gi,...,(?^;A,L) = [] Unh 

i<j ^ 

which coincides with that of the Chern-Simons model ()3.9p at A = 1 and without the Gaussian 
factors in the wavefunction. Then one obtains a relation analogous to (|3.5p , with an additional 
phase factor due to the appearance of hyperbolic sine functions instead of trigonometric ones 
(this factor also arises when one uses (|3.10p ). Repeating the procedure in ()3.2p using the relation 
sinh(ix) = i sinx we see that we now have to take L = —i (k + N) and that the additional phase 
factor, including the corresponding renormalization ()3.4p . is exp (^IttA^^), which also appears 
in the random matrix theory computation of ajsf |14j . 




3.3. Extension to generic root systems. The results we have described are not restricted to 
unitary gauge groups, and can be directly extended to orthogonal and symplectic gauge groups 
using generalizations of the Sutherland model appropriate to other semi-simple Lie algebras 
The generalized Calogero-Sutherland model is defined by the quantum Hamiltonian operator 



(3.16) ^ = -EI;i + E 



where A is the root system of a semi-simple Lie algebra in an A^-dimensional real vector space 
V with inner product (— , — ), q = [qi, . . . , qj\[) G V, and g\a\ = \a\ {^\a\ ~ 1) vr^/L^ are coupling 
constants which depend only on the lengths |a| of the root vectors a G A (i.e. on the orbits of 
the Weyl group). The usual Sutherland model (|1.9p corresponds to the ^at-i root system. 
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The BCn root system contains three distinct Weyl group orbits and hence three couphng 
constants. The Hamiltonian is 

1=1 ^« i<j ^ 



(3-17) + ^2 + 12 — 

i=l 2=1 

with corresponding ground state wavefunction 

N 

(3.18) Qn; \i, X2, X,L) =11 sm^-[^) ^^^'^ 

i=l 

X sm — sm ■ 



2L \ 2L 



The case Ai = corresponds to the Sat root lattice, A2 = to the Cn root lattice, and 
Ai = A2 = to the -Dat root lattice. This generalized Sutherland model is related, as in the 
case of the ^at-i system, to the partition function of Chern-Simons gauge theory on with 
orthogonal and symplectic gauge groups. The generic framing dependence is described by ()3.14p . 

For illustration, let us consider explicitly the Lie group Sp{2N). The configuration is again 
fixed and equispaced as 

(3.19) {qi,...,qN) = {N,N -l,N -2,...,1) , 

and this choice corresponds to the Weyl vector p of the Weyl chamber of C^r. Note that there is 
no essential difference with the ^tv-i system treated earlier, where the equispaced configuration 
was taken to be = c — i with arbitrary but fixed c G M. Evaluating (|3.18|) with A = Ai = ^ 
and A2 = yields the probability density [36] 



(3.20) (vI/o(iV,iV-l,...,l;i,0,i,L))'= J] 

where 



i V 



(3.21) /(j) 



iV - I - i j odd < 
iV - i j even < 
TV - i + i j odd > iV 
N - ^ + 1 j even > N . 

3 



With the identification ()3.3p . this is the partition function of Chern-Simons theory on S with 
gauge group Sp{2N) [55j. 

3.4. Derivation of Wilson line observables. The Sutherland model also delivers other quan- 
tum topological invariants, besides partition functions, in a similar way. For example, quantum 
dimensions [56j can be obtained in this way as well. In Chern-Simons gauge theory, the Wilson 
line invariants associated to the unknot give rise to the quantum dimensions which reduce, in 
the semiclassical limit A; — )• 00 (5(5 — )• 0), to the dimensions of representations of the gauge group. 
For example, the irreducible representations R of the gauge group SU {N) can be parametrized 
by the lengths of the rows of Young tableaux pi, i = 1, . . . , N , with ^1 > ^2 > • • • > Mat. Using 
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again (|3.4p and (|3.2p . but now with these weights as the positions qi = Hi of the particles, i.e. 
an inhomogeneous configuration, one obtains 

(3.22) dinig R = Tyjj(unknot) = J] sin ~ 

again with the identification (j3.3p . In this way knot and hnk invariants in Chern-Simons gauge 
theory correspond to non-equihbrium configurations of particle positions in the Sutherland 
model. It would be interesting to study excited states of the Sutherland model in its classi- 
cal equilibrium configuration to see if they give rise to quantum dimensions or other observables 
of Chern-Simons theory. 

3.5. Hamiltonian analysis. Recall that the hyperbolic Sutherland model and the Chern- 
Simons model differ in a Gaussian factor at the level of the ground state wavefunction p.9p . In 
the relationship with Chern-Simons theory, this is only a (framing) phase contribution to the 
partition function. At the level of the corresponding Hamiltonians, they are related by [13 

(3.23) Hcs = ^hyp + ^ E - '^i) (^^T^) 

where -ffhyp is the Hamiltonian of the hyperbolic Sutherland model. The reason for the appear- 
ance of a two-body term at the Hamiltonian level to explain a one-body factor in the ground 
state wavefunction is explained in [47j (together with some other properties of the model). 

At large distances, the two-body potential is given by a one-dimensional Coulomb potential 
V{x — x' ) = \x — x' \. Thus it has strongly confining properties, and while both the Sutherland 
model and the Chern-Simons model are connected to Chern-Simons gauge theory, the addition 
of such a potential leads to a very different correspondence. For example, in j57j the equivalent 
relationship between the Sutherland model and the Thirring model is presented, and a gener- 
alized Sutherland model with a two-body potential is also studied. It is shown there that in 
order for the correspondence between the Sutherland and Luttinger-Thirring models to hold, 
the generic two-body potential V{x) cannot decay slower than \x\~'^ for large x^ As explained 
in ^57j, the one-dimensional Coulomb potential has Fourier transform of the form V{k) = A/k"^ 
and the physics is that of one-dimensional quantum electrodynamics, rather than the Luttinger- 
Thirring behaviour of the Sutherland model. Indeed, a linear confining potential (such as the 
one-dimensional Coulomb potential) is one of the main ingredients of a massive Schwinger model. 
It would be interesting to understand further the implications of the fact that the two-body po- 
tential of the Chern-Simons model is, at large separations, a one-dimensional Coulomb potential. 
The exactly solvable Chern-Simons model is also related to other systems in condensed matter 
physics [33], such as Laughlin and multilayer wavefunctions on thin cylinders, and also through 
a connection between the Sutherland model and Luttinger liquids with boundaries. 
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